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ABSTRACT 

This  paper  is  an  application  of  the  Piecewise-Sinusoidal 
Reaction  Matching  Technique  (PSRMT)  to  linear,  parallel 
thin-wire  dipoles.   The  "method  of  moments"  if  briefly 
outlined  and  the  fundamental  electromagnetic  principles 
of  reciprocity,  surface  equivalence  and  reaction  are  dis- 
cussed.  Through  a  logical  combination  of  these  concepts, 
PSRMT  is  developed  and  applied  to  thin-wire  dipoles  radiating 
in  free  space  due  to  a  delta  function  voltage  generator  at 
their  center.   The  driving  point  impedance  and  the  current 
distribution  are  calculated  and  compared  to  results  obtained 
from  other  independent  theoretical  techniques. 

The  major  contribution  of  this  investigation  is  not  that 
of  presenting  a  new  technique  for  calculating  the  input  im- 
pedance and  current  distributions  on  dipole  antennas.   The 
study  is  primarily  concerned  with  comparing  other  well-known 
techniques  with  that  of  PSRMT.   Emphasis  is  placed  on  the 
ease  and  convenience  of  expressing  the  mutual  impedance 
integrals  of  coplanar  dipoles  in  closed  form. 
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I.   INTRODUCTION 

A.   DESCRIPTION  OP  ANTENNA  PROBLEM  UNDER  INVESTIGATION 

This  paper  is  an  investigation  of  the  efficiency  of 
the  Piecewise-Sinusoidal  Reaction  Matching  Technique  (PSRMT). 
In  the  solution  of  problems  using  the  "method  of  moments" 
[3],  computation  time  on  the  computer  depends  primarily 
upon  the  amount  of  computational  effort  required  to  generate 
the  elements  in  the  generalized  impedance  matrix.   Usually 
this  relates  to  the  amount  of  numerical  integration  required 
to  calculate  the  Z   's.   However,  by  employing  piecewise- 
sinusoidal  basis  functions  to  obtain  the  impedance  matrix, 
closed  form  expressions  can  be  obtained  for  the  integrations. 
Also,  the  size  of  the  impedance  matrix  required  to  adequately 
solve  the  antenna  problem  is  relatively  small  when  compared 
to  the  use  of  other  basis  and  weighting  functions.   For 
these  reasons  PSRMT  will  be  applied  and  the  results  compared 
to  other  techniques. 

The  antenna  configuration  under  investigation  will  be 
restricted  to  a  linear,  thin-wire  dipole.   The  antenna  is 
assumed  to  be  radiating  in  free  space  and  is  driven  at  the 
center  by  a  delta  function  voltage  generator.   In  this 
particular  application  of  PSRMT,  all  segments  of  the  fila- 
mentary dipoles  are  parallel  and  of  equal  length.   This  is 
not  a  necessary  condition,  but  one  of  convenience.   The 


driving  point  impedance  and  the  current  distributions 
will  be  calculated  by  an  application  of  PSRMT. 

B.   BRIEF  HISTORY  OF  RELATED  WORK 

In  the  past  much  work  has  been  done  in  the  area  of 
determining  the  mutual  impedance  between  two  elements  with 
an  assumed  sinusoidal  current  distribution.   Carter  [5] 
has  presented  the  mutual  impedance  between  two  wires  in 
echelon  where  each  antenna  is  an  odd  number  of  half-wave- 
lengths long.   The  mutual  impedance  expression  between 
two  identical  nonstaggered  parallel  center-fed  antennas 
of  Brown  [1]  is  of  historical  renown.    The  work  of  Brown 
was  carried  one  step  further  by  Cox  [6]  who  developed  the 
mutual  impedance  between  parallel  antennas  of  unequal  lengths 

In  1957,  H.  E.  King  published  a  paper  [10]  that  developed 
an  expression  for  the  mutual  impedance  between  two  staggered 
parallel,  center-fed,  infinitely  thin  antennas  of  unequal 
lengths.   Lewin  [12]  analyzed  the  mutual  impedance  of  two 
coplanar  half-wave  dipoles,  in  arbitrary  positions.   Finally 
Lewin' s  work  was  generalized  by  Richmond  and  Geary  [10] 
for  any  length  antenna. 

All  of  the  above  studies  were  restricted  to  assumed 
sinusoidal  dipole  current  distributions.   Without  assuming 
a  current  distribution  upon  an  antenna  the  current  distri- 
bution may  be  found  by  breaking  the  antenna  up  into  several 
segments,  calculating  the  mutual  impedance  between  each 
pair  of  elements  in  turn  and  generating  an  impedance  matrix 
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that  characterizes  the  current  distribution.   This  calcula- 
tion of  the  mutual  impedance  is  done  In  general  by  the 
expression  published  by  Richmond  and  Geary.   The  technique 
crudely  described  is  that  of  reaction  matching. 

C.   SCOPE  AND  LIMITATIONS  OP  THE  STUDY 

This  paper  develops  the  multisegment  reaction  technique 
from  fundamental  electromagnetic  theory.   Although  the 
technique  is  applied  only  to  linear  thin-wire  dipole's  of 
perfect  conductivity,  the  procedure  may  be  extended  to 
more  general  and  intricate  geometries.   To  cite  one  of 
the  many  examples,  PSRMT  has  been  successfully  employed 
by  Travieso-Diaz  [20]  to  the  analysis  of  antennas  and 
scatterers  with  arbitrary  surface  shape. 


11 


II.   OUTLINE  OF  THE  METHOD  OF  MOMENTS 

The  following  discussion  portrays  the  major  steps 
employed  in  the  method  of  moments.   The  object  is  to 
reduce  the  functional  equations  of  field  theory  to  matrix 
equations  which  can  be  readily  solved.   The  name  "method 
of  moments"  has  been  given  to  the  mathematical  procedure 
for  obtaining  the  matrix  equations  [8], 

A.  FORMULATION  OF  THE  OPERATOR  EQUATION 

Consider  only  equations  of  the  inhomogenous  type 
L(f)  =  g,  where  L  is  a  given  operator,  g  is  the  excitation 
or  source  (known  function),  and  f  is  the  field  or  response 
(unknown  to  be  determined). 

One  must  determine  from  physical  considerations  that 
the  problem  in  question  is  deterministic,  that  is,  if  the 
solution  for  f  when  L  and  g  are  given  is  unique.   The  method 
of  moments  gives  a  general  procedure  for  treating  field 
problems,  but  the  details  of  solution  vary  widely  with 
the  particular  problem.   The  analysis  to  follow  is  that 
which  is  germane  to  the  antenna  problem  under  investigation. 

B.  POSTULATION  OF  THE  BASIS  FUNCTION 
Consider  the  deterministic  equation 

LCD  -  g  ,  (i) 
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where  L  is  a  linear  operator,  g  is  known,  and  f  is  to  be 
determined.   Let  f  be  expanded  in  a  series  of  functions  or 
a  linear  combination  of  a  set  of  N  vectors  |f  V 


f  =  Ho      f    ,  (2) 

*—   n  n   ' 
n 


where  the  c  '  s  are  constants  that  are  not  known  but  are 
n 

to  be  determined  in  the  solution  process.   The  f  's  are 
chosen  to  be  linearly  independent  and  are  called  expansion 
functions  or  basis  functions.   For  exact  solutions,  Equation 
(2)  is  usually  an  infinite  summation  and  f  '  s  form  a  complete 
set  of  basis  functions.   For  approximate  solutions,  Equation 
(2)  is  usually  a  finite  summation.   Substituting  Equation 
(2)  into  Equation  (1),  and  using  the  linearity  property 
of  L, 

Z  cn  L(f  )  =  g   .  (3) 

n 

C.   POSTULATION  OF  THE  WEIGHTING  FUNCTION 

It  is  assumed  at  this  point  that  a  suitable  inner 
product  <^f,g>  has  been  determined  for  the  problem.   The 
inner  product  is  a  scalar  defined  to  satisfy 

<f\g>  ■  <&,& 
<*f  +  £g,h>  =  «<f,h>  +  ^g,h> 

<f*,f>  >  0,    if  f  /  0 

=  0,    if  f  =  0  , 
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where  tf  and  8  are  scalars,  and  *  denotes  complex  conjugate 
[8].   The  appropriate  inner  product  for  this  particular 
antenna  problem  will  be  developed  in  the  section  of  this 
paper  entitled  "Reaction  Concept".   Nov/  define  a  set  of  M 
linearly  independent  weighting  or  testing  vectors  \vi   j,  and 
take  the  inner  product  of  Equation  (3)  with  each  w   .   The 
result  is 


Zc   /w  ,Lf  >  =  <w  ,g> 


w 


m=  1,  2,  3,  •••  M,  where  the  N  c  's  are  unknown  and  all 
other  terms  are  known  [3]. 


D.   SOLUTION  OF  THE  LINEAR  ALGEBRAIC  SYSTEM 

Equation  (4)  can  be  written  in  matrix  form  as 


w  h  =  m 


(5) 


where 


IH 


w1,Lf1    wlfLf2 


•      • 


in 


N 


and 


h 


w1,g 


wOJS 


w3,g 


A  consistent  equation  results  by  requiring  M  and  N  to 

be  equal.   That  is,  the  basis  set  f  and  the  testing  set 

w   each  has  N  elements.   The  coefficient  matrix  is  deter- 
m 

mined  by  inversion: 


N 


I 


mn 


-1  r 


ni 


(6) 


provided  that  the  matrix 


I 


mn 


is  nonsingular.   The 


solution  for  f  is  given  by  Equation  (2): 


• 

■ 

f 

n 

c 
n 

a 

. 

(7) 


This  solution  may  be  exact  or  approximate,  depending  upon 

the  choice  of  f   and  w   and  the  extent  of  the  summation. 

n      m 

The  particular  choice  w   =  f   is  known  as  Galerkin's  method. 

m    n 

E.   COMMENTARY  ON  THE  METHOD  OF  MOMENTS 

One  of  the  most  critical  tasks  in  any  particular  problem 
is  the  choice  of  the  basis  and  testing  functions,  f  and 
w  .   The  f  ' s  should  be  linearly  independent  and  chosen  so 
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that  some  superposition  (see  Equation  (2))  can  approximate 

f  reasonably  well.   The  w  's  should  also  be  linearly 
°  m  ° 

independent  and  chosen  so  that  the  products  <^w  ,g^ 

depend  on  relatively  independent  properties  of  g.   Obviously 

additional  considerations  affecting  the  choice  of  f   and 

&  n 

w  are  a.)   the  accuracy  of  solution  desired,  b.)  the  ease 
m  ' 

of  evaluation  of  the  matrix  elements,  c.)  the  size  of  the 

matrix  that  can  be  inverted,  and  d.)   the  realization  of 

a  well-conditioned  matrix  IT     1[8], 

1  mnJ 

P.   PIECEWISE-SINUSOIDAL  BASIS  SETS 

The  physical  problem  may  be  represented  by  many 

different  operator  equations,  and  a  suitable  one  must  be 

chosen.   There  are  also  an  infinite  number  of  sets  of 

expansion  functions  f  and  testing  functions  w  that  may 
f  n  &  m        j 

be  chosen. 

Intuitively,  the  more  nearly  the  expansion  set  resembles 
the  actual  form  of  the  current  distribution,  the  fewer  the 
terms  required  to  accurately  represent  the  current.   In 
addition,  certain  types  of  basis  functions  may  lessen  many 
time-consuming,  numerical  operations  incorporated  in 
approximations.   The  piecewise-sinusoidal  basis  set  exhibits 
these  qualities.   The  objective  is  to  represent  the  current 
distribution,  over  individual  subsections  of  a  domain  with 
expansion  functions  that  possess  greater  curve-fitting 
potential  than  do  flat  pulses  and  yet  lead  to  well-conditioned 
matrices  in  approximate  solution  procedures.   The  application 
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of  piecewise-sinusoids  in  Galerkin's  method  and  the 
method  of  moments  to  electromagnetic  field  problems  is 
known  as  the  reaction  concept  [3,17]. 
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III.   THEORETICAL  MOTIVATION 

The  following  discussion  will  develop  the  basic  concept 
utilized  in  the  formulation  of  the  reaction  matching 
technique  of  J.  H.  Richmond.   The  reaction  concept  of 
Rumsey  [17] >  the  reciprocity  principle  [9]  and  the  surface 
equivalence  theorem  [18]  were  presented  by  Richards  [13]  and 
Travieso-Diaz  [20].   The  concepts  of  reciprocity,  reaction 
and  surface  equivalence  form  the  theoretical  motivation 
for  PSRMT;  therefore  they  are  included  again  here.   With 
a  logical  combination  of  these  elements,  one  can  develop 
a  general  formulation  of  electromagnetic  boundary-value 
problems,  in  particular,  reaction  matching. 

A.   RECIPROCITY  THEOREMS 

Consider  a  pair  of  time-harmonic  sources  of  the  same 
frequency  w   radiating  in  a  linear,  isotropic,  but  not 
necessarily  homogenous  medium  with  permittivity  t  and 
permeability  J*.   One  of  the  sources,  designated  as  source 
1,  has  impressed  electric  and  magnetic  current  densities 
(J,,  M  )  which  radiate  fields  (E, ,  H , )  in  the  presence  of 
a  second  source,  source  2,  acting  as  a  receiver.   Similarly, 
the  current  densities  (J2,  M ? )  of  source  2  set  up  fields 
(Ep,  Hp)  in  the  presence  of  source  1  acting  as  a  receiver. 

The  fields  of  these  sources  at  every  point  in  the 
medium  are  related  by  the  point  reciprocity  theorem  of 
Lorentz  [9]  as  follows: 

18 


V   •  (E1  X  H2  -  E2  X  H1)   =   (J1  •  E2  -  J2  •  Ex) 


-   (M1  •  H2  -  M2  •  H1) 

(8) 


Now  let  the  sources  be  of  finite  extent,  and  assume  all 
inhomogeneities  in  the  medium  to  have  a  finite  distribution. 
Define  non-overlapping  surfaces  S,  and  Sp  enclosing  the 
source  distributions  (J-,,  M,  )  and  (Jp,  Mp)  respectively; 
all  inhomogeneities  are  contained  in  the  volumes  V,  and  Vp 
bounded  by  S,  and  S?.   Therefore,  the  medium  outside  S, 
and  Sp  is  homogeneous.   Under  these  conditions,  the 
reciprocity  theorem  of  Carson  [4]  establishes  a  relation 
between  the  source  distributions    (J,  ,  M,  )  and  (Jp,  IYL) 
and  the  fields  (E  ,  H  )  and  (Ep,  H „ )  of  Equation  (8)  as 


1  2 


(9) 


B.   SURFACE  EQUIVALENCE  PRINCIPLE 

The  surface  equivalence  principle  can  be  briefly 
described  as  follows  [13].   First,  an  original  source 
distribution  (J. ,  M. )  set  up  fields  (E,  H)  everywhere  in 
some  defined  medium.   Next,  the  original  source  is  turned 
off,  but  simultaneously  surface  current  densities  (J  >  M  ) 

S      o 

are  set  up  on  a  surface  S  which  encloses  the  original 

19 


source  and  any  inhomogeneities .   If  the  surface  current 

densities  ( J  ,  M  )  are  related  to  the  original  fields 
s    s 

(E,  H)  by 


J   =  n  X  H  (10A) 


n     =  E  X  n    ,  (10B) 


A 

where  n  is  the  unit  normal  at  S  directed  outward  from  S, 
then  these  surface  current  distributions  will  radiate  the 
original  fields  (E,  H)  in  the  region  exterior  to  the  surface 
S  and  will  generate  a  null  field  inside  S. 

C.   REACTION  CONCEPT 

In  195^  V.  H.  Rumsey  published  a  paper  [17]  in  which 
he  introduced  a  new  physical  observable,  "reaction". 
Define  a  complex  number  denoted  by  <a,  b>   as  follows: 

<a,  b>   =     JjjCJCa)  •  E(b)  -  M(a)  •  H(b)]dv, 

(11) 

where  the  volume  V  contains  the  source  a. 
The  reciprocity  theorem  states  that 

<a,  b>   =  <t>,  a>  (12) 

provided  that  all  media  are  isotropic  and  that  a  and  b  can 
be  contained  in  a  finite  volume.   The  scalar  <^a,  b>  is 
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a  measure  of  the  reaction  (or  coupling)  between  the  sources 
a  and  b.   Equation  (11)  is  to  be  understood  as  a  formula 
for  the  measure  of  reaction.   It  is  obvious  that  Equation  (9) 
may  be  restated  as  Equation  (12),  where  the  inner  product 
represents  the  complex  integrals  of  Equation  (9).   Equation 
(12)  is,  therefore,  the  reaction  form  of  reciprocity. 
These  theorems  will  now  be  used  to  develop  the  "Piecewise 
Sinusoidal  Reaction  Matching  Technique"  (PSRMT)  of  J.  H. 
Richmond. 
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IV.   PIECEWISE-SINUSOIDAL  REACTION  MATCHING  TECHNIQUE  (PSRMT) 

A.  INTRODUCTION  TO  PSRMT 

In  this  section,  the  method  of  moments  and  the  reaction 
concept  will  be  applied  to  the  solution  of  the  current 
distribution  on  linear  antennas.   This  particular  class 
of  problem  involves  only  antennas  which  consist  of  or  can 
be  represented  by  a  collection  of  thin-wire  segments. 

An  axial  filamentary  test  source  will  be  developed. 
After  specifying  a  suitable  current  distribution  on  both 
the  original  and  the  test  antenna,  the  reaction  concept 
will  be  applied.   Reaction  will  be  employed  in  order  to 
develop  a  system  of  linear-  equations  whose  solution  will 
yield  the  complex  coefficients  that  complete  the  description 
of  the  approximate  current  distribution  specified. 

The  reaction  concept  and  its  applications  have  been 
discussed  by  Rumsey  [17]  and  Harrington  [9].   The  interpre- 
tation of  PSRMT  to  follow  is  primarily  that  of  Richards  [13] 
and  Leonard  Tsai  [21].   Most  of  the  material  presented 
here  has  been  previously  presented  by  Richmond  [1*1,  15,  16] 
and  Richards;   it  is  presented  here  for  completeness. 

B.  THIN-WIRE  PROBLEM 

Consider  the  thin-wire  dipole  shown  in  Figure  1.   The 
dipole  arms  are  wire  segments  of  length  L.  and  L„ ,  not 
necessarily  equal,  but  shorter  than  *>/2 .   The  radius  "a" 
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FI GURE  1 .     A  1 
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of  the  wire  is  arbitrary,  but  in  the  thin-wire  case  always 

small  in  terms  of  A;   also,  a  «  L   and  a  «  L? .   The 

dipole  is  driven  at  the  junction  of  its  arms  with  a  magnetic 

ring  current  M, . 
to  d 

Let  the  dipole  depicted  in  Figure  1  radiate  in  free 
space,  setting  up  time-harmonic  field  intensities  (E  ,  H , ) 
everywhere.   By  the  surface  equivalence  principle  this 
antenna  may  be  replaced  by  a  set  of  surface  currents 
(J,,  M , )  enclosing  the  primary  source  where 


Jd  =   n  X  Hd  (13A) 


Md   =   Ed  X  n  (13B) 


at  the  surface,  and  n  is  the  unit  normal  vector  pointing 
outward  from  the  surface.   Let  this  surface  exist  at  an 
infinitesimal  distance  outwardly  away  from  the  actual 
dipole,  so  that  the  dipole  itself  is  replaced  by  surface 
currents  (J,,  M-),   Then  it  is  known  that  exterior  to  this 
surface  the  two  fields  (E, ,  H,)  exist  and  there  is  a  null 
field  inside  the  surface.   Now  the  original  source  can  be 
removed  without  disturbing  the  fields,  so  that  there  is 
free  space  everywhere. 

A  test  source  is  placed  at  the  axis  of  the  actual 
dipole  in  Figure  1,  in  order  to  employ  the  reaction  test. 
Alone,  this  filamentary  test  dipole  radiates  the  fields 
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(E, ,  H . )   in  free  space.   The  Carson  reciprocity  theorem  of 
Equation  (9)  may  be  applied  between  the  test  dipole  and 
the  unknown  equivalent  surface  current  on  the  original 
dipole  to  give: 

JhJa   '   Et  "  Md  '  Ht]dse  =  J[It  •  Ed  -  Mt  •  Hd^  . 

e  (1M) 


where  S   is  the  surface  of  the  original  dipole  as  shown  in 
Figure  1,  I   is  the  coordinate  along  the  filamentary  test 
dipole,  and  (I  ,  M, )  represents  the  currents  of  the  test 
dipole.   Due  to  the  null  field  generated  by  the  original 
source  in  the  interior  of  S  ,  where  the  coordinate  along 
the  filamentary  test  dipole  is,  the  right-hand  side  of 
Equation  (1*0  is  identically  zero.   Therefore  the  reaction 
between  the  two  sources  is  zero, 


<a,  b>  S  <b,  a>  =  0   , 

hence  the  term  "zero  reaction  condition".   This  zero 
reaction  theorem  was  developed  by  Rumsey  [17].   Equation  (lM) 
written  in  light  of  the  zero  reaction  condition  is: 


U- 


Md  •  Ht  dse  "   sJJJd  •  Et  dse 

e  (15) 


For  a  typical  antenna  problem,  the  excitation  consists 
of  a  voltage  generator  driving  the  dipole  at  the  junction 
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of  its  arms.   Since  the  antenna  is  thin,  a  "delta  gap" 
voltage  generator  is  assumed.   Here  the  impressed  current 
densities  are: 


Jd  =   0  (16A) 


Md  =  0    (  £-  a,  I   -  lQ)Vd      ,  (16B) 


where  V,  is  the  strength  of  the  voltage  generator,  "a" 
is  the  radius,  £fi  is  the  value  of  the  coordinate  I   at  the 
gap  and  b   is  the  Dirac  delta  function. 

In  this  case  M,  is  a  ring  of  magnetic  current  at  the 
terminals,  and  the  integral  on  the  left-hand  side  of 
Equation  (15)  reduces  to  a  line  integral.   Applying  Ampere's 
law  to  perform  the  integration  around  the  gap: 


.  jK  •  Ht  dse   "   "a  Vd   J  Va'V 
'e 


r2rr 

S 


■  "Vd  h   (£o}  •        (17) 


Combining  Equations  (15)  and  (17)  yields 

qjJJd  '  Et  dse   "   -Vd  \    ^  (18) 

be 

where  I,  (-^q)  is  the  complex  terminal  current  on  the  test 
dipole. 
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For  the  evaluation  of  Equation  (18)  allow  the  thin-wire 
approximation  (a  «  "M  to  hold,  thus  the  distribution  of 
current  density  may  be  assumed  to  be  uniform  around  the 
cylinder.   This  assumption  may  be  expressed  mathematically 
as  simply,  J,  =  1,/2-rra  where  "a"  is  the  radius  and  I, 
varies  only  along  the  length  of  the  dipole.   Thus  Equation 
(18)  becomes 

e 

Equation  (19)  is  the  result  of  replacing  the  surface 

current  density  J,  with  an  equivalent  filamentary  current 

of  strength  I,  =  27raJ,   which  is  located  on  the  previous 

surface  where  J,  resides.   The  integration  6.1     is  along 

d  °         e 

this  line  on  the  dipole  surface  and  E.  represents  the  field 
of  the  test  dipole  on  this  line.   It  should  be  noted  that 
if  the  test  dipole  were  not  axially  located,  then  this 
simplification  would  not  be  judicious. 

Equation  (19)  forms  the  foundation  upon  which  the  two 
segment  reaction  matching  technique  is  based.   Once  V, 
is  specified  and  the  current  distribution  I,  is  assumed, 
the  coefficients  for  the  current  I,  may  be  readily  determined, 
For  the  remainder  of  this  paper  the  discussion  will  involve 
current  distributions  on  both  the  test  dipole  and  on  the 
original  dipole  that  are  "piecewise-sinusoidal"  in  form. 
This  will  lead  to  considerable  simplification  because 
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a.)   closed  form  expressions  for  the  E,  term  can  be  found, 
and  b.  )   closed  form  expressions  for   /  I -  •  E,  6.1        also 
exist  for  coplanar  dipoles.   The  basis  for  these  claims 
will  be  presented  after  a  generalization  is  made  for  N 
segment  reaction  matching. 

An  essential  feature  of  any  acceptable  test  source 
(J, ,  M . )  must  be  that  the  fields  it  radiates  be  known  and, 
if  possible,  that  they  be  known  in  closed  form.   One  such 
source  is  a  filamentary  electrical  source  with  time-harmonic 
current,  as  shown  in  Figure  2. 

Let  the  source  be  located  on  the  z  axis  with  endpoints 
z..  and  Zp,  and  the  current  distribution  be  of  the  form: 


I(z)   =   A  e 


Jkz   4-   B  e 


-jk: 


(20) 


where  A  and  B  are  complex  constants,  k  =  2ir/x,   and  the 

1  vt 
time  dependence   e^    is  understood  and  suppressed. 


C.   ELECTRIC  FIELDS  OF  A  PIECEWISE-SINUSOIDAL  FILAMENTARY 
SOURCE 

As  expression  for  the  z  component  of  the  free-space 

electric  field  of  the  filamentary  source  of  Figure  2  is 

given  by  Schelkunoff  and  Friis  [19]  as 


1 
TpTrJuTZ 


j  t    e _±    _   j  ,    e . 

1        R.  2        R^ 


Vv 

>z 


(21) 


He 


re    I1   =   I(z1),    Ij   =    I'(z1),    I2   =    I(z2),    IJ   =    I'(z2), 
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Figure  2.     Filamentary  electric  line  source  on  the  z  axis 
sets  up  a  field  (Ef  ,E  )   at  an  observation 
point  U,z). 
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and  V  is  the  electric  potential  of  the  end  charges;  £ 

is  the  permittivity  of  the  medium.   Throughout  this  paper 

the  prime  denotes  taking  the  derivative  with  respect  to  z. 

The  electric  field  E  is  composed  of  both  a  normal 
component  E q  and  a  tangential  component  E  as  depicted 
in  Figure  2. 

Note  that  in  performing  the  dot  product  "J,  •  E  " 
of  Equation  (18)  for  linear  antennas,  only  the  z  component 
of  the  electric  field  makes  a  contribution.   That  is  the 
justification  for  ignoring  the  E»  component  when  considering 
only  linear  segment  configurations. 

The  term  ^v/^z  will  be  dropped  immediately  since  it 
vanishes  for  adjacent  segments  and  the  contributions  from 
the  gradient  of  the  potential  V  due  to  end  charges  for  the 
end  segments  may  be  assumed  negligible  [13].   Note  that 
although  this  expression  for  E   is  relatively  simple,  this 
closed  form  expression  is  rigorous  even  in  the  near-zone. 

The  complex  coefficients  A  and  B  in  the  expression 
for  I(z)  (see  Equation  20)  can  be  solved  in  terms  of  the 
end  currents  I,  and  I?.   The  resulting  equation  for  I(z) 
is 


I(z)   =   [I.,  sin  k(z2  -  z)   +   I~  sin  k(z  -  z..  )  ]  /sin  kd 

(22) 


where  d  as  indicated  in  Figure  3  is  the  length  of  the 
segment . 
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FIGURE  3.  Piecewise-sinusoidal  current  distribution 
on  an  antenna  segment. 
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Figure  3  illustrates  the  current  distribution  expressed 
in  Equation  (22).   It  is  clear  that  from  Equation  (22) 
I'  and  Ip  may  also  be  expressed  in  terms  of  I.,  and  I?. 
Since 


I'(z)   =   [-kl,  cos  k(z2  -  z)   +  kip  cos  k(z  -  z^)]  /sinkd 

(23) 


it  follows  that 


1^   =   k(I2  -  I1   cos  kd)  /  sin  kd  (2*JA) 


I£   =   k(I2  cos  kd  -  I  )  /  sin  kd   .  (2*»B) 


Note  that  from  Equations  (21)  and  (24)  the  electric  field 
intensity  of  the  sinusoidal  electric  line  source  may  be 
expressed  entirely  in  terms  of  the  geometry  of  the  observation 
point  and  the  end  currents  I,  and  Ip. 

D.   EXPANSION  OF  THE  CURRENT  DISTRIBUTION 

The  forms  of  the  current  distributions  on  both  of  the 
filamentary  dipoles  is  piecewise-sinusoidal.   The  current 
on  each  segment  of  each  dipole  can  be  described  by  the 
sinusoidal  functions  given  in  Equation  (22).   Because  the 
components  of  the  electric  field  intensity  of  such  an 
element  are  known,  both  terms  within  the  integral  of  Equation 
(19)  exist  as  closed  form  expressions. 
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The  following  piecewise-sinusoidal  distribution 
has  been  employed  by  Richmond  [15]: 

-,  i         tfcl   A         sin  k(l-l      _  ) 

j   =    1  lU)      =   1    t  1  *   P  (£)  n- 

s     2-rta  2ta   ^in  n  n  n 


n=l       1X        sin  kd 


a  sin  k(£  .,  -  I)      /oc.s 

+  T  4.-,  P  .i.-,  W — (25) 

n+1   n+1  sin  kd 


or 


s     2tia 


N-l 

I 

n=l 


=   ~-   7   I   G  (26) 

Tta   *-,   n  n 


with  I   =  1(1    )    and  d  the  length  of  the  segment.   P  (I) 
n      n  to  to        n 

is  a  unit  pulse,  nonzero  only  on  segment  n;  f   is  a  unit 
vector  pointing  in  the  direction  of  I   from  one  endpoint  of 
segment  n  to  the  other.   The  vector  functions  G   are 
known  as  basis  functions,  expansion  functions  or  dipole 
modes.   The  mode  currents  go  to  zero  at  the  endpoints 
of  the  dipoles,  as  depicted  in  Figure  H. 

It  is  advantageous  to  have  a  symmetric  impedance  matrix 

Z    (to  be  defined  later).   Furthermore,  the  test  sources 

mn  *■ 

should  be  chosen  to  yield  a  well-conditioned  system  of 

linear  equations  for  the  subsequent  inversion  process. 

For  these  reasons  and  to  obtain  a  closed  form  expression 

for  the  integral  of  Equation  (19),  the  expansion  function 

for  the  test  source  G  will  be  of  the  same  form  as  the 

m 

basis  function  for  the  original  dipole' s  current  distribution 
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FIGURE  4.  Basis  functions  (Gn)  representing  the 
assumed  current  distributions  on  a 
linear  filamentary  antenna. 
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With  the  piecewise-sinusoidal  current  distribution 
employed  for  each  segment  of  the  antenna,  then  from 
Equation  (22)  the  current  on  any  given  segment  can  be 
considered  the  sum  of  two  sinusoidal  currents,  each  of 
which  vanishes  at  one  or  the  other  end  of  the  segment. 
This  may  be  visualized  with  the  aid  of 'Figure  H.      One  can 
see  how  the  equivalent  surface  distribution  is  treated  as 
a  collection  of  (N-l)  overlapping  tubular  dipole  current 
densities  on  an  antenna  consisting  of  N  segments.   The 
nth  dipole  has  terminal  current  I   at  the  nth  junction  with 

r  n  o 

the  current  vanishing  at  the  far  ends  of  the  two  segments 
forming  the  nth  junction. 

The  solution  for  the  unknown  complex  terminal  currents 
Si  \    proceeds  by  introducing  a  set  of  (N-l)  filamentary 
dipoles  each  having  piecewise-sinusoidal  currents  and 
being  located  along  the  axis  of  the  original  antenna.   No 
generality  is  lost  by  assuming  each  test  dipole  has  unit 
terminal  current.   Reciprocity  is  enforced  successively 
between  each  test  dipole  and  all  the  tubular  surface 
dipoles  in  turn.   The  "zero  reaction  condition"  yields 
the  following  equation  for  the  mth  test  dipole  and  the 
nth  surface-current  dipole  (refer  to  Equation  (15)): 


(27) 


The  left-hand  side  of  Equation  (27)  is  nonzero  only 
when  n  =  m.   This  is  the  result  of  the  tangential  component 
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n 

M      ■ 

•    H      ds 
m        n 

— 

n 

|Jn' 

•    E      ds 
m        n 

of  the  electric  field,  on  the  surface  of  the  original  source 

being  zero  except  at  the  delta  gap.   M  which  equals 

E   x  n  is  therefore  also  zero  unless  the  mth  test  dipole 
n  K 

is  co-located  with  the  nth  original  dipole.   Expressed 
mathematically  : 


/M   •  H  ds   =   -V     /  H  dl     £  =   -V  I  £ 

n    mn       n   »J  m   n  mn       nm 

n 


nj     iu   n     um  n  m  mn 

(28) 


The  right-hand  side  of  Equation  (27)  can  be  simplified 
by  carrying  out  the  circumferential  integration  of  J   •  E  . 
Again  this  is  equivalent  to  replacing  the  tubular  distri- 
bution with  filamentary  segments  for  calculating  reaction: 


|f  J   •  E   ds    =    /  I  (I    )    •  E  dl 

JJ  n    m   n      »J  n   n     m   n 


(29) 


n  n 


Equating  Equations    (28)    and    (29)    one   obtains 


-V     T     i  =  [i    (i    )    •    E      d/         .  (30) 

nmmn  »    J     nMn  mn 


/„ 


v  <r 

= 

- 

I    (T   )-E      dT 

m     mn 

vn 

n      n        m        n 

Assuming  unit  amplitude  for  the  test  current  and  noting 
that  upon  co-locating  m  =  n, 

(3D 

The  integration  is  carried  out  on  t  ,  the  coordinate  along 

to  n»  o 

the  arms  of  the  nth  surface  dipole  filaments.   The 
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symbol  $        is  the  standard  Kronecker  delta,  defined  by 
J      mn  s  j 

the  relations: 


i  =0    m  J   n 

mn 


yi     0    m  =  n 


(32) 


Equation  (31)  forms  the  basis  of  the  multi-segment  reaction 
matching  technique. 

E.   MATRIX  EQUATION  FORMULATION 

Substitute  the  assumed  form  of  I  (i-  )  (refer  to  Figure 

n  n  to 

H)   namely 


I  (r  ) 
n  n 


=  I 


n 


sin  k(£  -  I     _ ) 


sin  kd 


sin  k(£  .,  -  £) 

+  *  xi  P  A1U)  -^ 

n+1  n+1        sin  kd 


(33) 


into  Equation  (31)  and  define  a  mutual  impedance  term  as 


mn 


-  -  I 

n-1 


sin  k(£  -  £   , ) 

*^-T   •  E  d^ 

sin  kd      n    m 


-./ 


n+1   sink(£n+1  -  I)   . 


n 


sin  kd 


rn+l-Em  «■  (3<4) 


The  dummy  variable  £  is  the  coordinate  along  the  particular 
segment,  the  other  variables  are  defined  as  in  Equation  (25). 
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E   as  before,  is  the  free  space  electric  field  intensity 
m         ' 

radiated  by  the  mth  test  dipole  with  piecewise-sinusoidal 
current,  specifically  for  the  mth  test  dipole  composed  of 
a  lower  segment  ML  and  an  upper  segment  MU,  the  z  component 
of  the  electric  field  is 


-jkR 


1ML 


1ML 


.  J'kR2ML 
'I2ML  e 


m 


H-njv*t 


R 


1ML 


R 


■JkR 


1MU 


2ML 


•jkR 


2MU 


1MU 


2MU 


R 


1MU 


where  as  depicted  in  Figure  5,   !•! 


R 


2MU 


(35) 


1ML 


=   I 


'   of  segment 


ML,  I-Jjvttt   =   I-,  of  segment  MU,   rimt   =   Rn  of>  segment  ML, 


R 


1MU 


R..  of  segment  MU   etc... 


Figure  5  is  merely  an  extension  of  the  single  segment 
with  piecewise-sinusoidal  currents  of  Figure  2  to  a  dipole 
composed  of  two  such  segments. 

Equation  (31)  may  now  be  written 


V  i 

m  mn 


mn  n 


(36) 


where  I   is  the  complex  terminal  current  at  the  nth  junc- 

n  * 

tion.   The  calculation  of  the  reaction  between  the  mth 
test  dipole  and  the  entire  equivalent  surface  current 
density,  consisting  of  N-l  overlapping  surface  dipoles, 
yields  the  following  relation: 
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FIGURE  5.  The  m   test  dipole,  composed  of  a  lower  segment 
ML  and  an  upper  segment  Ml) 
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N-l 

V   =  T    z   I   .  (37) 

m     /TZj   mn   n  J 


Equation  (37)  represents  one  linear  equation  out  of  a 
system  of  (N  -  1)  equations  that  are  obtained  by  successively 
enforcing  reciprocity  (or  calculating  reaction)  between  the 
(N  -  1)  surface  dipoles  and  each  of  the  (N  -  1)  axial  test 
dipoles  in  turn.   The  test  is  performed  in  turn  between 
each  of  the  axial  dipoles  and  the  tubular  dipoles  yielding 
a  system  of  linear  equations  that  can  be  described  by  the 
matrix  relation 


[Z   ]  [I  ]   =   [V  ]   ,  (38) 

mnJ   n        m  J   * 


which  is  of  the  same  form  as  Equation  (37),  except  that 

Z  _  is  now  a  square  (N  -  1)  X  (N  -  1)  matrix  and  V  is  a 
mn  ^  m 

column  vector.   This  set  of  equations  can  be  solved  for 
the  unknown  junction  current  coefficients  \1   \     after  a 
known  voltage  source  is  placed  at  the  actual  feed  junctions 
and  the  voltage  of  the  other  junctions  are  set  to  zero. 

The  coefficients  I   of  Equation  (38)  can  be  solved  for 
by  several  techniques.   Crout  published  a  paper  [7]  in 
which  he  discusses  a  method  for  evaluating  determinants 
and  solving  systems  of  linear  equations  with  complex 
coefficients.   In  the  case  of  well  conditioned  matrices 
the  impedance  matrix  can  be  merely  inverted  and  multiplied 


ho 


by  the  voltage  matrix  to  yield  the  terminal  current 
coefficients : 

-1 

[I  ]   =   [Z   ]    [V  ]  (39) 

n        mn      m  VJ^y 

Equation  (39)  may  be  written  simply  as ' 


-1 
[I]   =   [Z]    [V]  (1I0A) 


or  equivalently : 


[I]   =   [Y]  [V]   ,  (l|0B) 


where 


-1 
[Y]   =   [Z] 

The  complex  current  coefficients  I   that  are  determined 
from  Equation  (38)  are  used  to  establish  the  current 
distribution  found  from  Equation  (26).   Once  the  current 
distribution  is  known,  parameters  of  interest  such  as  field 
patterns,  input  admittances,  etc.,  can  be  calculated  by 
numerically  evaluating  the  appropriate  formulas. 

A  wire  antenna  is  obtained  when  the  wire  is  excited 
by  a  voltage  source  at  one  or  more  points  along  its  length. 
For  an  antenna  excited  at  the  nth  junction,  the  applied 
voltage  matrix  of  Equation  (40)  is 


ill 


[V  ]   = 
n 


V 


n 


(4l) 


that  is,  all  elements  of  the  voltage  matrix  are  zero  except 
the  nth,  which  is  equal  to  the  source  voltage.   The  current 
distribution  is  established  by  Equation  (40)  which  for 
the  [V  ]  of  Equation  (4l)  becomes 


[I]   =   V. 


n 


- 

Yln 

Y2n 

Y3n 

« 

• 

• 

Y 
mn 

. 

(42) 


Note  that  the  nth  column  of  the  admittance  matrix  gives 
the  current  distribution  for  a  unit  voltage  source  applied 
to  the  nth  junction.   Inversion  of  the  impedance  matrix 
therefore  gives  simultaneously  the  current  distributions 
when  the  antenna  is  excited  at  any  arbitrary  junction 
along  its  length.   For  the  case  where  only  one  junction 
is  excited  with  a  voltage  V. ,  the  input  admittance  is 
given  directly  by  the  feed  junction  current  I.  as 
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YIN   "   vf   '  ^3) 


P.   SOLUTION  TECHNIQUES 

As  mentioned  previously,  integrals  of  the  form  of 
Equation  (3*0  can  be  obtained  in  closed  form  for  coplanar 
segments.   In  the  case  of  non-coplanar  segments,  the 
integrals  must  be  calculated  by  numerical  methods,  for 
example,  Simpson's  rule  or  Gaussian  quadrature. 

Upon  a  careful  examination  of  Equations  (3*0  and  (35) 
it  should  be  ovbious  that  the  integral  of  Equation  (35) 
may  be  expressed  in  closed  form  with  the  relationship  [16] 
that  follows: 


rv        exp(ijr)  exp(-jR  ) 

/  d      —  dr   =   [exp(-jz  )  E(R  -z  +r) 

J  „R  m     m  m 


m  vi  rR 

1  m 


r=r 

-  exp(jz  )  E(R  +z  ?r)]   2 
v       m     m  m    r=r. 


W) 


where 


Ci(x)   =  V  +  Ln(x)   +  J     cos  (t)  '  1  dt    (45A) 


Jsin(t 
t 


Si(x)   =    r   SiSill  dt   -  |    .  (*45B) 
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Utilizing  the  relationship  in  Equation  (44),  Richmond 
and  Geary  published  a  paper  [16]  that  developed  a  closed 
form  expression  for  the  mutual  impedance  between  coplanar 
skew  dipoles.   However,  for  this  antenna  problem,  all  of 
the  reactions  will  be  performed  between  parallel  dipoles. 
A  more  useful  closed  form  expression  for  the  problem 
consisting  of  entirely  parallel  dipoles  was  presented  by 
H.  E.  King  [10].   The  complete  expression  with  its  deriva- 
tion may  be  found  in  Appendix  A  and  B.   Even  though  these 
closed  form  expressions  are  not  particularly  simple,  they 
do  involve  only  sine  and  cosine  integrals  which  have 
polynomial  algorithms;  hence  Z   in  this  form  is  more 
efficient  than  "brute-force"  integration. 

Leonard  Tsai  [21]  formally  portrayed  PSRMT  as  a 
particular  application  of  Galerkin's  method.   Consider 
again  the  operator  equation 

L(f)   =   g  (46) 

where  f,  the  given  operator,  is  specifically  the  unknown 
current  distribution  on  a  filamentary  dipole  and  the 
operator  L  yields  the  fields  g  of  a  distance  "a"  away 
(on  the  original  dipole).   Expanding  f  into  a  linear 
combination  of  basis  functions: 

f   =  E  c   f  (47) 

n     ~*       n   n 
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and   enforcing  Galerkin's   method  directly: 


V     c        <Lf    ,    f    >     =    <g      f    >        .  (48) 

*-»       n  n*      irr  >      m  ' 

n  ' 


Here  the  f  would  describe  the  weighting  functions 
corresponding  to  the  "test  dipoles"  which  is  using  the 
same  expansion  functions  as  the  "original  dipoles",  both 
subdomain  piecewise-sinusoidal  pulses.   The  inner  product 
of  the  left-hand  side  of  Equation  (48)  corresponds  directly 
to  the  integral  of  Equation  (3D: 

I  •  E  dr  (49) 

+  J      n    m   n 

n 

used  in  the  reaction  matching  technique.   Note  also  that 

<g,  f  >  =  V  which  is  nonzero  only  at  voltage  sources. 
>OJ   mm  *         o 

This  follows  from  the  "zero  reaction  condition".  Thus 
the  procedures  of  PSRMT  are  those  employed  by  a  direct 
application  of  Galerkin's  method. 
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V.   DISCUSSION  OF  NUMERICAL  TECHNIQUES  AND  COMPUTER  PROGRAM 

Evaluation  of  the  integrals  in  Equation  (3*0  may  be 
accomplished  by  any  one  of  several  numerical  integration 
techniques  such  as  Romberg,  Weddle,  Gaussian  Quadrature, 
etc. ,  as  well  as  by  the  closed  form  expression  developed 
in  Appendix  A.   For  most  integrals,  no  representation  in 
terms  of  elementary  functions  is  possible,  and  approximation 
becomes  necessary.   Most  numerical  techniques  are  successfully 
based  upon  some  type  of  polynomial  approximation:   trapezoi- 
dal rule,  Simpson's  rule,  Romberg's  method,  etc..  Unfortu- 
nately, these  techniques  have  the  usual  error  sources  such 
as  truncation,  but  perhaps  more  interestingly  they  are 
harassed  with  the  question  of  convergence.   The  question 
is  whether  as  continually  higher  degree  polynomials  are 
used  or  as  continually  smaller  intervals  between  data 
points  are  used,  a  sequence  of  approximations  is  produced 
for  which  the  limit  of  the  truncation  error  is  zero.   As 
may  be  guessed,  the  smaller  the  interval  h  is  made,  the 
more  computation  time  required  and  the  greater  the  rounding 
off  error  becomes.   In  other  words,  accuracy  is  limited 
for  any  given  method  by  algorithm  errors  that  ultimately 
obscure  convergence. 

It  is  found  in  practice  that  decreasing  the  interval 
between  sampling  points  below  a  certain  level  leads  to 
larger  errors  rather  than  smaller.   This  necessitates 
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experimentation  of  the  various  parameters  available  in 
the  technique  in  order  to  establish  a  trend  (if  one  exists) 
that  gives  acceptable  results.   Unfortunately  this  requires 
an  assumed  correct  solution  in  order  to  evaluate  the 
performance  of  the  integration  on  a  supposedly  unknown 
integral. 

Table  1  illustrates  how  the  value  for  input  admittance 
of  a  half-wave  dipole  varies  as  a  function  of  the  number 
of  antenna  segments  and  the  order  of  the  Gaussian  Quadrature 
used  in  evaluating  the  integral  of  Equation  (3^0  •   Table  2 
is  a  similar  portrayal  of  a  full-wave  antenna. 

The  norm  for  Table  1  and  2  was  chosen  to  be  the  value  of 

input  admittance  given  by  King  [11].   Table  1  shows  no 

convergence  tendency  whatsoever,  while  Table  2  indicates 

that  a  lower  degree  polynomial  may  be  used  in  evaluating 

Z    if  the  antenna  is  divided  into  more  segments.   The 
mn  to 

apparent  discrepancy  between  Table  1  and  2  is  probably 
explained  by  the  difference  in  the  actual  current  distri- 
butions on  the  two  antennas.   It  is  clear  that  the  problem 
is  multi-faceted. 

This  convergence  problem  is  not  exclusive  to  the 
Gaussian  Quadrature  technique.   In  fact,  all  of  the 
numerical  techniques  have  some  troublesome  parameter  analysis 
involved  in  order  to  effectively  use  them.   But  a  far  more 
important  and  exciting  reason  to  avoid  their  use  is  that 
of  efficiency.   The  closed  form  expression  involves  only 
the  evaluation  of  the  sine  and  cosine  integrals  which  have 
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TABLE  1.  Convergence  analysis  for  the  input  admittance 
of  a  half-wave  dipole.  The  numbers  respresent 
percent  variation  from  the  King-Middleton  standard 
of  Y=  10.8  mmho's.  (Percentages  expressed  to  the 
nearest  .5%) 
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TABLE  2.     Convergence  analysis  for  the  input  admittance  of  a 
full -wave  dipole.     The  numbers  represent  percent 
variation  from  the  King-Middleton  standard  of 
Y  =  -90  mmho's.      (Percentages  less  than  20  expressed 
to  the  nearest  .5%) 
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simple  algorithms.   Unfortunately  the  closed  form  solution 

is  valid  in  general  for  skewed  coplanar  dipoles  only, 

therefore,  numerical  integration  would  be  required  .for 

non-coplanar  elements.   The  same  antenna  problem  (if 

coplanar)  may  be  solved  by  the  closed  form  technique  in 

a  fraction  of  the  time  that  numerical  methods  require. 

The  computer  program  listed  in  Appendix  B  calls  upon 

subroutine  CZMN2  for  the  evaluation  of  the  impedance 

matrix.   It  is  here  that  the  method  of  solution  may  be 

chosen.   If  the  closed-form  route  is  decided  upon, 

CZMN1  (which  is  also  listed)  is  used  to  evaluate  Z 

mn 
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VI.  ANALYSIS  OF  RESULTS 

A.   CONVERGENCE  CHECK 

Figures  6  and  7  illustrate  the  accuracy  and  efficiency 
of  PSRMT.   Leonard  Tsai  [21]  used  similar  figures  (the 
solution  labeled  "Richmond's  Reaction")  to  show  the  compara- 
tive rate  of  convergence  of  PSRMT  to  that  of  point  matching 
with  triangular  pulses.   Note  that  reaction  matching 
required  only  six  segments  while  for  the  same  accuracy 
point  matching  required  five  times  as  many  segments. 
Superimposed  on  Tsai's  comparison  are  those  points  generated 
by  the  closed  form  expressions  of  Equation  (72)  of  Appendix 
B. 

Figures  8  and  9  illustrate  the  efficiency  of  PSRMT. 
In  Figure  8,  as  more  and  more  flat  pulses  are  employed 
in  Galerkin's  method  the  more  the  solution  corresponds  to 
that  of  PSRMT  with  only  six  segments.   Figure  9  shows  how 
rapidly  PSRMT  relates  to  the  moment  method  employing 
piecewise-sinusoidal  basis  functions. 

Figure  10  exemplifies  the  convergence  of  PSRMT  itself. 
The  comparison  is  drawn  between  PSRMT  and  Hallen's  equation 
employing  piecewise-sinusoidal  basis  functions.   16  segments 
do  an  admirable  job  in  relating  to  this  half-wave  dipole. 
This  is  sufficiently  true  under  certain  conditions  (as 
will  be  shown)  to  warrant  a  rule  of  thumb:   use  32  segments 
per  wavelength. 
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Figure  11  is  a  plot  of  conductance  versus  susceptance 
for  thin  cylindrical  dipoles  with  radius  a./%   =  .001588. 
The  locus  spirals  in  as  the  length  of  the  dipole  increases. 
Note  that  for  small  lengths  the  8  segmented  solution  devi- 
ates from  the  accepted  King-Middleton  value  of  admittance. 
Many  of  the  plots  to  follow  will  use  King-Middleton  values 
for  comparison,  which  were  taken  from  tables  published  by 
R.  W.  P.  King  [11].   In  the  evaluation  of  the  admittance, 
the  King-Middleton  procedure  was  again  used  (two  iterations 
and  an  empirical  normalization  were  also  employed).   Figure 
11  indicates  the  advantage  of  taking  32  segments  per 
wavelength. 

Figure  12  is  the  same  type  as  Figure  11  but  the  radius 
has  been  increased  to  .00635^.   It  is  evident  that  for  short 
dipoles  with  large  radii,  the  length  to  radius  ratio  for 
each  segment  exceeds  the  thin-wire  approximation.   That  is 
why  that  as  the  length  of  the  antenna  is  increased  in  Figure 
12,  the  solution  becomes  more  convergent. 

Figures  13  and  14  portray  the  convergence  trends  of  the 
Weddle  numerical  integration  technique  on  a  half-wave  dipole 
As  indicated  in  Figure  13,  the  solution  varies  only  slightly 
with  a  change  in  the  number  of  integration  points  used  in 
the  integration.   However,  the  input  admittance  approaches 
the  King-Middleton  solution  most  nearly  when  the  number 
of  integration  points  equals  29.   Figure  14  shows  the  change 
in  current  distribution  as  the  number  of  antenna  segments 
change  (with  the  number  of  integration  points  equal  to  29). 
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Again  the  input  admittance  approaches  the  King-Middleton 
solution  as  the  number  of  segments  increase.   These 
particular  trends  are  valid  only  for  this  half-wave  dipole 
of  radius  a/\  =   .001588.   Other  geometries  must  be  considered 
on  a  one  to  one  basis  due  to  the  behavior  of  Weddle's 
technique  on  the  integral  of  Equation  (3*0. 
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FIGURE  6.   Input  resistance  versus  h/7.  for  a  cylindrical  wire 
antenna  with  radius  a  and  total  length  L  =  2h. 
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FIGURE  7.  Input  reactance  versus  h/-\  for  a  cylindrical  wire 
antenna  with  radius  a  and  total  length  L  =  2h. 
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Galerkin  Method  with  Flat  Pulses 


.  _  „  ,. —  _  5  pU]ses 

— — 13  pulses 

— 20  pulses 

30  pulses 

PSRMT  -  Closed  Form 

A       A       *      6  segments 


FIGURE  8.     Current  distribution  on  a  half-wave  dipole  with 
radius  a/x  =   .001588. 
(Galerkin  Method  with  flat  pulses) 
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General  Moment  Method 


6  basis  functions 

-13  basis  functions 

20  basis  functions 

30  basis  functions 


PSRMT  -  Closed  Form 


6  segments 


FIGURE  9.     Current  distribution  on  a  half-wave  dipole  with 
radius  h/x  =   .001588. 
(General  moment  method  -  piecewise-sinusoidal   pulses) 
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Hall  en's  Equation 
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Input  Impedance 
86.62  +  j  46.78    /v 


PSRMT  -  Closed  Form 
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32 

92.25  +  j   37.38    „_ 

FIGURE  10.     Current  distribution  on  a  half-wave  dipole  with 
radius  a/*  =   .001. 
(Hallen's  Equation  -  piecewise  -  sinusoidal   -  basis  functions) 
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King  -  Middleton 
PSRMT  -  Closed  Form 


^   8     segments 
o      32L/>.  segments 
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(mmho's) 


(mmho's) 


FIGURE  11.     Admittance  curve  for  thin-cylindrical   di poles  with 
radius  a/%  =    .001588. 
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PSRMT  -  Closed  Form 
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+  16  segments 
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X  32  segments 
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FIGURE  12.  Admittance  curve  for  thin  cylindrical  dipoles 
with  radius  a/x  =  .00635. 


60 


Symbol 
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PSRMT  No.  of  Pts 

(numerical    integration)   5 
Neddie's  Method  17 
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FIGURE  13.     Current  distribution  on  a  half-wave  dipole  with 
radius  a/v  =   .001588  (convergence  check  on  the 
number  of  integration  points  used  in 
Weddle's  technique) 
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FIGURE  14.     Current  distribution  on  a  half-wave  dipole  with  radius 
a/x  =   .001588. 

(convergence  check  on  the  number  of  segments  used  in 
Weddle's   technique) 
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B.   ADMITTANCE  CORRELATION 

Figures  15  through  18  attempt  to  typify  the  admittance 
correlation  between  the  King-Middleton  solution  and  various 
PSRMT  solution  techniques.   Figure  15  and  16  depict  the 
solution  by  Weddle's  method  and  the  Gaussian  Quadrature 
method,  respectively.   Figures  17  and  18  can  be  compared 
with  Figure  11  to  show  the  effect  of  changing  the  dipole's 
radius.   The  effect  is  negligible  as  long  as  the  length 
to  radius  ratio  of  each  segment  is  maintained  large. 
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King-Middleton  » 

PSRMT  -  Numerical    Integration       a 
(Weddle's  Method 
Number  of  Integration 
points  equals  29) 


(mmho's) 


FIGURE  15.     Admittance  correlation  for  thin  cylindrical   dipoles 
with  radius  a/%  =   .001588.    (32  L/*  segments) 
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6th'  order) 
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FIGURE  16.  Admittance  correlation  for  thin  cylindrical  dipoles 
with  radius  a/x  =  .001588.   (32  L/x  segments) 
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FIGURE  17.     Admittance  correlation  for  thin  cylindrical   dipoles  with 
radius  a/*  =   .003175.      (32  L/tc segments) 
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FIGURE  19.  Admittance  correlation  for  thin  cylindrical  dipoles 
with  radius  a/x  =  .004763.  (32  L/>  segments) 
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C.   CURRENT  DISTRIBUTIONS 

Figures  19  through  25  present  current  distributions  on 
dipole  antennas  of  various  lengths  and  radii.   On  each  set 
of  curves  is  indicated  the  expansion  set  used  to  represent 
the  unknown  current.   The  number  of  segments  employed  in 
the  PSRMT  technique  was  based  upon  the  rule  of  thumb 
32  LA  . 

In  Figures  24  and  25,  the  power  series  [2]  used  as  a 
comparison  may  be  expressed  as 


£  c  zn- 


Kz)   = 

ri=l 


where  the  C  's  are  complex  constants  to  be  determined. 
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King-Mi ddleton 


PSRMT  -  Closed  Form        A 
(8  segments) 


FIGURE  19.     Current  distribution  on  a  quarter-wave  dipole  with 
radius  a/%  =   .001589. 
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King-Mi ddleton 


PSRMT  -  Closed  Form         A 
(8  segments) 


FIGURE  20.     Current  distribution  on  a  quarter-wave  dipole  with 
radius  a/>  =   .007022. 
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King-Mi ddleton 


PSRMT  -  Closed  Form  a 

(16  segments) 


FIGURE  21.     Current  Distribution  on  a  3/4-wavelength 
dipole  with  radius  a/%  =   .001588. 
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King-Middleton 


PSRMT  -  Closed  Form 
(32  segments) 


FIGURE  22.     Current  distribution  on  a  full -wave  dipole 
of  radius  a/*  =   .001583. 
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King-Mi ddleton 


PSRMT  -  Closed  form 
(48  segments) 


FIGURE  23.     Current  distribution  on  a  3/2    "Mipole  of  radius 
a/>  =   .001588. 
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King-Mi ddleton 

Butler-Power  series 

PSRMT  -  Closed  form 
(16  segments) 
/>  -  Gaussian  Quadrature 

(4  segments) 


FIGURE  24.     Current  distribution  on  a  half-wave  dipole  with 
radius  a/%  =   .001588. 
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FIGURE  25.     Current  distribution  on  a  half-wave  dipole  with 
radius  a/%  =   .005. 
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VII.   CONCLUSIONS  AND  RECOMMENDATIONS 

A.   EVALUATION  OP  POTENTIAL  INTEGRALS 
Integrals  of  the  form 


rh         -jRl       "JR2 
F(z)   =      f(z')  [S .   +  £ ]  dz»         (50) 

<y  r1  r2 


where 


R1  =  "^/(z-z')2  +  a2 


R2  =   V(z+2')2  +  a2 


are  extremely  difficult  to  evaluate  for  "small  a"  because 
of  the  sharp  peak  in  the  integrand  at  z'  =  z.   Although 
PSRMT  has  been  formulated  on  the  thin-wire  assumption 
(a  <<"*-);   for  very  small  radius  values  the  integrand  in 
Equation  (3*0  becomes  ill-behaved  due  to  the  form  of  the 
electric  field  E   (refer  to  Equation  (35)). 

This  phenomena  was  observed  in  the  calculation  of  the 
input  impedance  of  Figures  6  and  7.   In  fact,  for  h/x  <  .3 
the  value  of  the  radius  (h/a  =  2000)  became  small  enough 
to  thwart  computation  entirely. 
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B.  APPLICATION  OF  THE  GENERAL  PSRMT  TECHNIQUE 
Although  this  paper  dealt  in  particular  with  parallel 

dipoles  only,  PSRMT  is  useful  in  solving  a  wide  variety 
of  thin-wire  antenna  problems.   With  the  aid  of  the  concepts 
of  reciprocity,  reaction  and  surface  equivalence,  this 
technique  reduces  an  integral  equation  to  a  set  of  simul- 
taneous linear  equations.   This  set  of  equations  is  readily 
solved  on  a  computer  to  obtain  the  input  impedance  and 
current  distribution  on  the  antenna. 

Richmond  and  Geary  [16]  presented  a  new  form  for  the 
rigorous  near-zone  field  of  the  linear  dipole  with  sinusoidal 
current  distribution.   This  electric  field  expression  was 
employed  to  derive  a  closed  form  expression  for  the  mutual 
impedance  between  coplanar-skew  dipoles.   By  employing 
their  expression  for  mutual  impedance,  complex  geometries 
and  surfaces  may  be  represented  by  wire  grids  and  reaction 
enforced  to  yield  solutions. 

C.  RECOMMENDATIONS 

The  theory  developed  in  this  paper  is  by  no  means 
restrictive  to  the  specific  problem  investigated.   The 
data  comparisons  were  made  by  employing  PSRMT  to  equally 
segmented,  thin-wire  linear,  parallel  dipoles.   Extensions 
to  this  problem  are  obvious.   First,  the  segments  need  not 
be  of  equal  length.   By  merely  altering  the  computer  program, 
the  specific  segment  lengths  can  be  applied  to  the  appro- 
priate formulas  to  generate  a  more  general  solution.   Second, 
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the  dipoles  need  not  be  parallel.   By  employing  Richmond 
and  Geary's  closed  form  expression,  the  mutual  impedance 
between  any  two  coplanar  elements  may  be  determined. 
But  now  the  reaction  between  the  filamentary  test  dipole 
and  the  equivalent  filamentary  tubular  surface  dipole  must 
be  examined  in  more  detail.   In  the  case  of  parallel 
elements,  the  reaction  applied  between  the  test  dipole  and 
any  filamentary  surface  dipole  is  the  same.   But  due  to 
the  geometry  of  the  problem  involving  two  non-parallel 
segments,  the  mutual  impedance  depends  upon  the  angular 
position  of  the  filamentary  surface  dipole  with  respect 
to  the  axis.   The  question  of  how  does  one  best  react  these 
segments  is  yet  to  be  answered.   E.  G.  Neely  will  present 
a  Master's  Degree  Thesis  to  the  Electrical  Engineering 
Department  at  the  U.  S.  Naval  Postgraduate  School  (in 
September  1973) j  in  which  he  deals  with  arbitrarily  bent 
monopoles  above  a  ground  plane.   He  is  approaching  this 
problem  by  averaging  the  reactions  taken  around  the 
circumference  of  the  tubular  dipole  successively.   Other 
approaches  might  be  investigated. 
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APPENDIX  A 

TRANSFORMATION  OF  EXPONENTIAL  INTEGRALS 
INTO  CLOSED  FORM  EXPRESSIONS 

The  z  component  of  the  electric  field  from  a  dipole 
with  a  time-harmonic  current  distribution  on  each  segment 
is  expressed  by  Equation  (35)-   For  the  mth  test  dipole, 
the  dipole  mode  current  distribution  is  shown  in  Figure  26 


79 


!iMr° 
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FIGURE  26.  Typical  mode  current  distribution  for  a 
filamentary  test  dipole. 
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From  the  expressions  for  the  derivatives  of  the  current 
at  the  endpoints  of  the  segments  (given  in  Equation  (24)), 
one  obtains 


I1ML   =   k(I2ML  "  I1ML  COS  kd)  f    Sin  kd       (51A) 


X2ML   =   k(I2ML  C0S  kd  "  I1ML)  7  Sin  kd       (51B) 


I1MU   =   k(I2MU  "  I1MU  C0S  kd)  /  Sin  kd       (51C) 


I2MU   =   k(I2MU  C0S  kd  "  I1MU)  /  Sin  kd       (51D) 


Note  that  from  Figure  26,  it  is  apparent  that 


I1ML  =   0  (52A) 


X2ML  =   I1MU  =  X  C52B) 


I2MU  =  0  (52C) 


The  result  of  substituting  Equation  (52)  into  Equation  (51) 
is 


I1ML   =   -J2MU   =  k  f   Sln  kd  (53A) 


I2ML   =   "I1MU   =   k  COS  kd  f   Sin  kd  (53B) 
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Redefine  the  distance  designations  of  Figure  5  according 
to  Figure  27. 
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FIGURE  27.     Geometry  of  two  parallel   dipoles   in  echelon 
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With  the  geometry  indicated  in  Figure  27  and  Equation 
(53)  the  electric  field  of  Equation  (35)  may  be  rewritten 


as 


z 


-JkR. 


m 


4itjwt  < 


k   e 


k  e 


-JkR. 


R1    sin   kd 


-   2k   cos   kd 


Rp    sin   kd 


-JkR, 


Rn    sin   kd 


(5^) 


or  equivalently : 


E~     -      J  ClJiwI 


m 


k        s     A 

)    z 


-JkR. 


k  e 


k  e 


-JkR, 


R..    sin   kd  R~    sin  kd 


-JkR, 


+      2k   cos   kd 


Rfi    sin   kd 


(55) 


In  this  problem 


k  =  2-ir  /  *     with   "X  =   v/f  =   1 
a/  =  2-r(3(108)) 
£  =  1/36 -ir  (10"9) 


therefore 


k/%w*  =   30 


(56) 
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Equation  (55)  in  light  of  Equation  (56)  is 


E 


m 


j    30  | 

iin   kd 


<  — 


-JkR. 


R. 


-\jkR2  ~J'kRn 

e  e 

-  +    2    cos   kd 


R, 


R 


0 


(57) 


where 


■v: 


RQ       a2  +  z2 


(58A) 


i  -v 


a2  +  U  -  z)2 


R, 


Va2  +  ^i  +  z)2 


(58B) 


(58C) 


Employing  the  trigonometric  identity 


sin  2  °<  =   2  sin  «  cos  <»c   , 


Equation  (57)  may  be  written  as 


E 


m 


30  z 
sin  kd 


-J  e 


-JkR. 


R. 


-JkR2              -JkRQ 
-j  e 2j  cos  kd  e 


R, 


R 


0 


(59) 


Equation  (59)  is  exactly  the  form  that  H.  E.  King  expressed 
his  electric  field  in;   except  for  the  extra  sin  kd  term  in 
the  denominator.   This  term  is  not  present  in  King's 
expression  for  the  electric  field  because  he  was  assuming 
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that  the  antenna  current  distribution  is  sinusoidal. 
In  the  application  of  PSRMT  the  current  distributions  are 
assumed  to  be  piecewise-sinusoidal  and  therefore  the 
expression  for  the  mutual  impedance  is: 


mn 


-30 


(sin  kd) 


-J  e 


-JkR. 


k  U   -  I     ^ 


f     J       n  sin 
L  C  *n-l 

r*n+i  T 

+   J  sin  *Un+rD    f 


R. 


-jkR2 

,   -J  e ,   2j  cos  kd  e 

R2  RQ 


-JkR 


Equation  (60)  follows  from  the  substitution  of  Equation  (59) 
into  Equation  (3*0  and  assuming  that  the  length  of  all 
segments  are  equal  to  d. 

In  order  to  facilitate  transforming  Equation  (60)  into 

a  closed  form  expression.  Z   will  be  expressed  in  terms 

v  3      mn  ^ 

of  the  geometry  of  Figure  27.   For  sinusoidal  current 
distributions, 


Z    =   -30 
mn 


£0+h 


2l^+h 


♦      /""2 


sin  k(z  -  h) 


l2+h 


sin  k(2£2+h-z) 


-J  e 


-jkR. 


R, 


-JkR.  -JkR  \ 

-J_e_  ±   +   2J  cos  k  i  d 


R, 


R, 


dz 
61) 
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Obviously  Equation  (61)  can  be  multiplied  out  to  yield 
six  expressions,  one  of  which  follows: 


=   +30    J 


2l2+h  i  -jkF^ 


sin  k(2£„  +  h  -  z)   +l-£- 


l2+l  2  1     Rl 


dz 


(62) 


Employing  the  trigonometric  relations: 


\CX- 


eJ    =   cos  oc  +  j  sin  p<-  (63A) 


sin^   =  -sin  (-*)  (63B) 


cos  *   =  cos  (-oc)  (63C) 


Equation  (62)  becomes 

/ 


£2    /j  cos  kR    sin  kR 

a  =   30   J  5 -  +  — = -I  sin  k(2£0+h-z)dz 

l2+h  \        Rl         Rl 

(61]) 


The  real  part  of  complex  Equation  (6l)  gives  the  mutual 
resistance  and  the  imaginary  part  gives  the  mutual  reactance 
The  real  part  of  Equation  (6*0  gives  only  one-sixth  of  the 
total  expression  for  the  mutual  resistance,  but  it  will  now 
be  transformed  into  a  closed  form  expression  in  order  to 
illustrate  the  procedure  of  obtaining  all  twelve  terms, 
both  resistive  and  reactive. 
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/    2    /sin  kFL  \ 

AREAL   =   3°     *  — R Sln  k(2^2  +  h  "  z)  dZ 

(65) 


Substitute  the  expression  for  R_  given  in  Equation  (58) 


to  obtain 


r   2£_+h   ,   .->  /  2  _,_  ,       x2 
/    2   sin  k  ya  +  (    -  z) 


A  =  30   J        ■  sin  k(2£0+h-z)  dz 

REAL     U+h  -,  I   2  "  /       ^  2 


'2'11       "^a2  +  (  ~  -  z) 


(66) 


By  a  change  in  the  variable  of  x  =  I     -   z ,  Equation 
(66)   becomes 


£..-2£„-h 


ARPfiT    =   -30     J  si_n_kVA     +   x        sin  kC2^9+h-£.+x)    dx. 

REAL        Wh       I  VTT^     /  2     x 

(67) 


Employing  the  trigonometric  identity  for  the  product  of 
angles  one  obtains 
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/    1    2  -\T? 2 

+15     J  [-cos  k(h-£..+2£9)  cos  k(  a  +x   -x 


EAL     '^ih 


-  sin  k(h  -  £,  +  2^2)  sin  k(  a2+x2  -  x) 


+  cos  k(h  -  £  +  2l2)    cos  k(  a2+x2  +  x) 


sin  k(h  -  I.    +  2£0)  sin  k(V+x2  +  x)3  ■     dx 

1  d  V~T~2 

a  +x 


(68) 


A  second   change   of  variable   of 


u  =   k(a  +x2   -  x)  (69A) 


v  =   k(    a2+x2   +  x)  (69B) 


gives 


ru3 

AREAL   =    15   J  [  +  C0S    k(h    "    ll    +    2V    £27Ji      dU 

Ul 

+sin   k(h   -   I,    +    2^0)    sin   u]    du 
1  2  u 


f  3 

+   15      J         C  +  cos   k(h    -   £      +    2£2)    cojL-V      dv 

vl 


-sin  k(h   -  £1   +   2£2)    ^—^]    dv  (70) 
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where 


u  =  ^"^a  +  (h  -  I     +  l2)2   +  (h  -  £x  +  £2)) 


u3  =  k(%2  +  (h  -  l1   +    2l2)2   +    (h   -   l±   +   2^2)) 


vn  =  k(   a2  +  (h  -  £_  +  £0)2  -  (h  -  I     +  £Q)) 
1  12  i  2 


v3  =  k(   a2  +  (h  -  l±   +  2£2)2  -  (h  -  ^  +  2^2)) 


The  integrals  of  Equation  (70)  are  recognized  as  sine 
and  cosine  integrals  and  may  be  written  as 


AREAL  =  15   C0S  k(£l  "  2l2    "  h)  ^-Ci^u^  -  ci(vx)  +  Ci(u3)  + Ci(v3)] 


+  sin  k(£  -  2l2   -  h)  [Si(u1)  -  Si(v  )  -  Si(u  )  +  Si(v.)] 

(71) 


Recall  that  Equation  (6l)  through  Equation  (70)  was  written 
for  sinusoidal  current  distributions.   To  make  them  appli- 
cable to  piecewise-sinusoidal  current  distributions, 

2 

Equation  (6l)  must  be  multiplied  by  the  factor  1/  (sin  kd)  . 

Equation  (6l)  expanded  in  closed  form  as  typified  by 

4 

Equation  (71)  is  written  in  its  entirety  in  Equation  (72) 

and  (73)  in  Appendix  B.   Equation  (73)  is  an  expression 

for  the  mutual  impedance  between  two  thin  parallel  center-fed 
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antennas  with  half-lengths  HL  spaced  a  distance  DL  apart 
and  staggered  by  the  value  SL.   This  expression  is  a 
simplification  of  H.  E.  King's  general  formula  for  mutual 
impedance  [10]. 
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APPENDIX  B 
FORTRAN  PROGRAM  LISTING 

Figure  28  is  the  main  computer  program  listing  (in 
Fortran  IV)  for  the  solution  of  the  linear  dipole  antenna's 
driving  point  admittance  and  current  distribution. 

The  symbols  for  the  input  data  are  defined  as: 

EL   length  of  the  dipole  in  wavelengths 

A    radius  of  the  dipole 

N    number  of  antenna  segments. 

This  program  has  built  in  "do-loop"  that  generates 
the  voltage  matrix  of  Equation  (4l).   The  voltage  source 
of  one  volt  is  automatically  placed  at  the  center  junction 
of  the  dipole. 

Calculation  of  the  complex  Z  „  matrix  is  done  in  the 

subroutine  CZMN2.   This  calculation  may  be  performed  by 

one  of  several  techniques.   Gaussian  Quadrature,  Romberg 

and  Weddle's  methods  can  be  employed  in  order  to  effect 

numerical  integration  of  Equation  (3*0  for  Z   .   If  this 
&  M  mn 

is  done,  additional  subroutines  are  required  to  calculate 
the  electric  field  explicitly,  establish  the  limits  of 
integration,  perform  the  necessary  integration  technique, 
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etc..  Z   may  be  most  efficiently  calculated  by  programming 

directly  the  closed  form  expression  of  Equation  (73). 

Once  Z   was  determined,  CMIN1,  a  subroutine  developed 
mn  '       '  * 

by  E.  G.  Neely,  was  used  to  calculate  the  complex  inverse 

of  Z   .   The  complex  current  coefficients  result  from  a 
mn 

direct  application  of  Equation  (40).   Since  a  single  one 
volt  source  was  placed  at  the  center  junction,  the  driving 
point  admittance  is  known  immediately  from  Equation  (^3). 

The  program's  output  is  set  up  to  provide  a.)   an  echo 
check  on  the  input  parameters  and  important  variables, 
b.)   list  the  real  and  imaginary  parts  of  the  complex 
junction  currents  as  a  function  of  z/h,  c.)   list  the 
magnitude  and  phase  of  junction  currents,  and  d. )   display 
the  input  admittance  and  impedance.   This  output  is 
presented  for  both  the  dipole  and  the  monopole. 
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